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We outline the basic questions that are being studied in the theory of entanglement. Following a brief review 
of some of the main achievements of entanglement theory for finite-dimensional quantum systems such as 
qubits, we will consider entanglement in infinite-dimensional systems. Asking for a theory of entanglement in 
such systems under experimentally feasible operations leads to the development of the theory of entanglement 
of Gaussian states. Results of this theory are presented and the tools that have been developed for it are applied 
to a number of problems. 
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INTRODUCTION 

When John Bell, then a theoretical physicist at CERN, pub- 
lished his now famous theorem in 1964, it was at first hardly 
noticed by the scientific community. Only within the sub- 
sequent decade the fundamental character of its content was 
truly appreciated 1 1]: in a sense, it fleshed out the notion of 
entangled quantum systems being potentially correlated in a 
way that classical systems cannot. Entanglement can indeed 
be viewed as the furthest and most radical departure of quan- 
tum mechanics from classical physics. In recent years, yet an- 
other shift in the way quantum correlations are perceived has 
taken place: entanglement has been recognized as a resource 
that may be exploited to perform completely novel informa- 
tion processing tasks or enhance the efficiency of known tasks. 
As for any resource, a detailed understanding of its properties 
and the resulting possibilities and limitations for its manipu- 
lation are an important pre-condition for the full exploitation 
of its potential. Quite naturally, the systematic investigation 
of quantum entanglement in particular is a major goal of the- 
oretical research in quantum information science QElBSl- 

This program of research has been initiated with the inves- 
tigation of the entanglement properties of bi-partite quantum 
systems prepared in pure quantum states 13, la. Today, most 
characteristics of p ure state entanglement are in fact well un- 
derstood | §1 0, ll(ifllll . However, our understanding of mixed 
state entanglement is much less complete. Most proven re- 
sults are restricted to situations where the constituent parts are 
qubits, quantum two-level systems, which is due to a whole 
range of special mathematical properties only satisfied by 
qubit systems 1 12||. In contrast, the study of bi-partite entan- 
glement between higher-dimensional systems is much more 
involved. Indeed, new types of entanglement emerge such 
as bound entanglement |13], and one encounters a rich and 
complex proliferation of types of entanglement with growing 
dimension of its constituent parts. 

Considering the complex structure that is emerging even in 
the case of low-dimensional quantum systems, one may be 
tempted to come to the conclusion that to say anything mean- 
ingful about the entanglement of infinite-dimensional quan- 
tum systems such as field modes of light or harmonic oscilla- 
tors is an enterprise doomed to failure. This is, however, not 



so: it has become clear in recent years that generic statements 
on the entanglement of such quantum systems can indeed be 
made. Most importantly, for a large class of states many ex- 
act results on entanglement theory can be found, even in cases 
where there is no finite-dimensional counterpart. A brief de- 
scription of the development of this theory will be the main 
subject of this article. It should be noted that in this article, a 
certain emphasis is put on the results that the research groups 
at Imperial College and in Potsdam have contributed to this 
area Q. 

This article will be structured in the following way. After 
a short introduction to some key notions of entanglement that 
are valid for any dimensions we will then move on to outline 
the key questions that any theory of entanglement will aim 
to answer when it addresses the resource character of quan- 
tum entanglement. Some major results and definitions from 
the theory of entanglement for finite dimensional systems will 
be presented. From there we will move on to prepare the 
ground for the study of the entanglement properties of infinite- 
dimensional systems by introducing some basic notions that 
are helpful for the description of states and dynamics. Then 
we will exhibit some of the problems that occur when one 
attempts to study entanglement of infinite-dimensional sys- 
tems without any restrictions. This will serve as a motivation 
to restrict entanglement theory to so-called Gaussian states. 
Again, basic results are presented; we will then ask the exper- 
imentally motivated question of the manipulation of Gaussian 
states by Gaussian operations. We will outline results that 
have been obtained on this question, including a no-go the- 
orem concerning entanglement distillation of Gaussian states 
and Gaussian operations, and results on general entanglement 
manipulation for Gaussian states. In the conclusions we will 
outline further developments and applications of the theory of 
Gaussian entanglement in other areas of quantum information 
science. 



BASIC NOTIONS OF ENTANGLEMENT 

The power of entanglement manifests itself particularly 
clearly in quantum communication tasks that have to be per- 
formed over large distances. Such tasks will typically require 
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that, say, two parties involved in the communication protocol 
establish a shared entangled state of two quantum particles be- 
tween them. This may for example be accomplished by one 
party preparing two entangled particles locally followed by 
the transmission of one of these particles through a physical 
channel such as an optical fiber to the other party. For large 
distances this task is complicated by the unavoidable presence 
of unwanted interactions with the environment causing deco- 
herence and dissipation that will tend to destroy entanglement. 
The two distant parties that wish to create a shared entangled 
state between them will therefore generally end up with a par- 
tially entangled mixed state. Consequently, in order to be able 
to implement an quantum communication protocol, they will 
need to be able to 'repair' the partially entanglement state. 
This aim will be hampered by the fact that the parties are spa- 
tially separated, as this implies that they can only implement 
general quantum operations in their respective local laborato- 
ries and coordinate their respective actions by classical com- 
munication. Global quantum operations affecting particles in 
separate laboratories are not available to them. The set of local 
operations and classical communication is usually abbreviated 
LOCC and forms the basis for much of the study of quantum 
entanglement, though more general classes of operations can 
also be relevant 11 51 1 1 611 . The very notion of entanglement is 
actually tied to the set of LOCC: a state is called disentangled 
if it can be prepared using LOCC only. Indeed, the most gen- 
eral two-party state that can be generated from a product state 
employing LOCC is a mixture of product states of the form 

P = {Pi ® P? ) i (!) 

i 

where the tensor product refers to the parties A and B. Such 
a state, usually called a separable state, may exhibit correla- 
tions, but they are not of genuinely quantum mechanical na- 
ture as, again, this state may be created using LOCC UJ- A 
separable state allows for a description in terms of a local clas- 
sical model. Any state that cannot be cast into the form Eq. 
(or appropriately approximated by such a state for infinite 
dimensional systems) will be called entangled. The develop- 
ment of a theory of entanglement, or in other words the study 
of what can and cannot be achieved under LOCC, is a key con- 
cern in the field. Such a theory will generally aim to answer 
the following three central questions: 



Characterize 

While the definition of separability is easy to formulate, it 
is very difficult to decide in practice whether a given state is 
separable or not. Following Eq. Q, in order to show that a 
state is separable, it appears that one has to construct explic- 
itly a decomposition of the state into tensor products. This is 
a very difficult and potentially lengthy task especially for high 
dimensional systems. For low dimensional systems, however, 
the separability question can be decided in a different and 



more efficient way using the theory of positive but not com- 
pletely positive maps. In fact, a simple necessary and suffi- 
cient criterion for the separability of a quantum state can be 
based on the properties of the transposition and its application 
on a single sub-system. Clearly, the transposition is a positive 
map in the sense that it maps any positive operator onto a posi- 
tive operator, i.e., if p is positive then so is p T . The same then 
applies when the transposition is applied to one sub-system, 
say system B, of a separable state, because 

i 

is again a valid state. However, when we apply this so-called 
partial transposition to an inseparable state, then there is no 
guarantee that the result is again a positive operator, i.e., a 
physical state. Indeed, when applying the partial transposi- 
tion to the state with state vector \ip) = (|01) - |10))/v / 2 
(with respect to the basis {|00>, |01>, |10>, |11»), we find 
= |(|01)(01| + |10)(10| - |00)(11| - |11)(00|), 
which is evidently not positive any more. This observation 
gave rise to the conjecture that a state is inseparable exactly 
if the partial transpose is not a positive operator. It has in- 
deed been proven that this conjecture is correct for systems 
composed of two qubits or a qubit and a qutrit 1 1 8 ] . It should 
be noted that the partial transposition on qubits corresponds 
to time reversal. This observation already suggests the ap- 
propriate application of the partial transposition criterion for 
infinite dimensional systems. In systems with larger dimen- 
sion that six the situation is fairly different, and the question 
of separability is far more complicated: one encounters for 
example inseparable states whose partial transposition is nev- 
ertheless positive again - a phenomenon commonly referred 
to as bound entanglement I13l Il9ll . While there is a system- 
atic strategy for deciding the separability of a state 12011 . it is 
known that the problem has a computational complexity that 
liesinNP-HARD|2l]]. 

Manipulate 

A natural next question would be: Given a quantum state 
of which one knows that it is entangled, how can it be ma- 
nipulated by LOCC? One may envision, for example, the sit- 
uation that an experimental procedure generates a particular 
state p, but that a state p' different from p is actually required. 
Is it possible to achieve the transformation p — > p' employ- 
ing LOCC only? In general, this is a extraordinarily difficult 
question to answer, but some results have been achieved for 
finite dimensional systems. Most progress has been made 
for pure quantum states for which the basic structures are 
known. There are two regimes which may fruitfully be con- 
sidered. Firstly, one may study the situation where a single 
copy of a quantum state is held by two parties (the multiple 
copy case is a special case of this setting). In this case, the 
mathematical structure underlying the manipulation of quan- 
tum states by LOCC is provided by the theory of majorization, 
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and both necessary and sufficient conditions for the intercon- 
version between two states are known, together with protocols 
that achieve the task lliHsSl Holl . Let us write the initial and final 
state vectors as \tpi) and \4>2) in their Schmidt-basis, 

n m 

Hi) = E ^ m i**) > w») = E i*a> i4> , (3) 

where n denotes the dimension of each of the quantum sys- 
tems. We can take the Schmidt coefficients to be given in 
decreasing order, i.e., a\ > oli > . . . > a n and a[ > a' 2 > 
... > a' n . The question of the interconvertability between 
the states can then be decided from the knowledge of the real 
Schmidt coefficients only. One finds that a LOCC transfor- 
mation that converts to \1jj2} with unit probability exists 
if and only if the {a*} are majorized l! 22ll by {o^}, that is, 
exactly if for all 1 < I < n 

1 1 

X><XX (4) 

i=l i=l 

Various refinements of this result have been found that provide 
the largest success probabilities for the interconversion be- 
tween two states by LOCC together with the optimal protocol 
where a deterministic interconversion is not possible fl fioll . 
These results allow in principle to decide any question con- 
cerning the LOCC interconversion of pure state by employing 
techniques from linear programming [10]. Although the basic 
mathematical structure is well understood, surprising and not 
yet fully understood effects such as entanglement catalysis are 
still possible lHHl 

Secondly, further to the study of the LOCC transformation 
of single copies, a fruitful area is the study of the asymptotic 
regime ie the limit of a very large number of identical quantum 
states. While this may be viewed as a limiting case of the 
above theorem, it has indeed predated it. It turns out that in 
this case, a single number, the entropy of entanglement, will 
determine what is possible and what is not. The entropy of 
entanglement for a pure state with state vector \if>) is defined 
as 

E(MM) = (Sott B )(MM). (5) 

where S denotes the von-Neumann entropy defined as S(p) = 
— tr[p log 2 p\ and tr^ denotes the partial trace over subsys- 
tem B. With this quantity, we find that in the asymptotic 
limit of a large number of identically prepared systems, the 
transformation from N copies of \ip{) (t/)i\ under LOCC to 
M = NE^^D/E^iH) copies of \ip 2 )(H will be 
possible with fidelity that approaches unity in the limit of N 
approaching infinity. 

Quantify 

For mixed states no comparable set of theorems is avail- 
able. Indeed, even very simple questions concerning the in- 
terconvertability of states cannot be provided with necessary 



and sufficient conditions for their availability. The quest for 
such conditions together with the questions of the best effi- 
ciencies for procedures of salient interest that are available 
has led to the development of the theory of the quantification 
of entanglement. An entanglement measure grasps this no- 
tion of quantifying entanglement: it is a mathematical quan- 
tity that captures the essential properties that we would as- 
sociate with the amount of entanglement and that is ideally 
related to some operational procedure. Any function E map- 
ping the state space on positive real numbers that deserves the 
name entanglement measure should satisfy the subsequent list 
of requirements: 1E3H24I1 

1. E(p) vanishes if the state p is separable. 

2. E is invariant under local unitary transformations, i.e., 
local basis changes. 

3. E does not increase on average under LOCC, i.e., 

E{p)>Y^PiE{pi) (6) 

i 

where in a LOCC the state with label i is obtained 
with probability pi. 

4. For pure state | ip) (ip | the measure reduces to the entropy 
of entanglement 

£(MM) = (Sott B )(MM). (7) 

A function E mapping state space on positive numbers that 
satisfies the first three conditions is called an entanglement 
monotone while an entanglement measure also satisfies con- 
dition (4). Often, convexity of E or asymptotic continuity are 
taken as additional reasonable properties (see, e.g., Ref. 1 24]). 
In fact, conditions (1-3) together with asymptotic continuity 
uniquely specify the entanglement measure to be the entropy 
of entanglement on pure states (4), up to a real number that is 
merely a scaling of the quantity [ 24] . 

These conditions alone, yet, do not uniquely specify a mea- 
sure of entanglement on mixed states. The extremal measures 
J25I1 givin g bounds for all others are the distillable entangle- 
ment II15LI26II . quantifying the degree to which maximally en- 
tangled states can be extracted under LOCC from a given sup- 
ply of identically prepared states, and the entanglement cost 
lt26ll27ll . quantifying the resources in terms of maximally en- 
tangled states necessary in a preparation procedure of a state 
under LOCC. The entanglement cost is in fact nothing but the 
asymptotic version of the entanglement of formation I26ll28ll . 
The relative entropy of entanglement quantifies to which ex- 
tent agiven state can be distinguished from a separable state 
1E3I I29L I30I . In its asymptotic version l3lll it provides more- 
over a tight upper bound for distillable entanglement. Prac- 
tically most straightforward to use is probably the negativity 
132 1 13 311 . which is known to be an entanglement monotone 
II34II35II . and whose logarithm (up to a constant) has an inter- 
pretation as an asymptotic entanglement cost [16]. The very 
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recently developed 'squashed' entanglement |36] is a convex 
entanglement monotone that is moreover additive, a property 
that it shares with the relative entropy of entanglement with 
reversed arguments 1 37], but it remains finite on pure states. 
The same set of questions namely, characterization, manipula- 
tion and quantification will have to be addressed in the devel- 
opment of a theory of entanglement. In the remainder of this 
article we will outline what has been achieved in this respect 
in the infinite dimensional setting. 



THE TRANSITION FROM FINITE TO INFINITE 
DIMENSIONS 

All the above results concern bi-partite entanglement for 
systems comprising finite dimensional sub-systems such as 
qubits or qutrits. More recently, however, considerable ef- 
fort has been directed towards the study of infinite dimen- 
sional quantum systems such as the photon number degree of 
freedom of light modes |38, 39, 40, 41] or nano-mechanical 
harmonic oscillators l42ll . or the state of cold atomic gases 
l43ll . At first sight, one would expect the theory of entangle- 
ment to become extraordinarily complicated due to the fact 
that now the Hilbert space is no longer finite-dimensional. In- 
deed, without imposing certain restrictions on the set of states 
under consideration, one even loses such elementary proper- 
ties as the continuity of entanglement and its measures |44] 
and questions such as distillability become difficult to answer 
I45I1 . In this subsection we will demonstrate some of the coun- 
terintuitive properties I44I l46l l47ll of entanglement measures 
in the infinite dimensional setting and use them to motivate 
possible restrictions on the set of states that one wishes to con- 
sider. 

To clarify this issue, let us again consider the joint system 
to be bi-partite, consisting of parts labeled A and B, each of 
which having a finite number of degrees of freedom. Both Tt a 
and Hb of the joint Hilbert space Ti — Ha <g) 7Y_b are taken 
to be infinite-dimensional, as in case of two modes of light. 
We will consider trace-norm continuity, where the trace norm 
||.||i is defined as IU||i = tr[|A|] = tr[(A+ A) 1 ' 2 ] for trace- 
class operators A 14811 . It is now quite straightforward to see 
that the entropy of entanglement Q is no longer trace-norm 
continuous. The following sequence of states exemplifies this 
property of the entropy of entanglement: the entropy of entan- 
glement may be very different from zero or even infinite for 
states that are arbitrarily close to pure product states. In this 

example, let er = IV'oKV'ol witri IV'o) = \4$) ® \<t>is)> and 
define {<Jk}kLi as a sequence of pure states = |V'fc)(V'fc| 
defined by 

m = (i- £fc ) 1/2 hM 



(8) 



and e k = l/log(fc) 2 . Here, {|</£°) = n £ N } and 
{ 1 4>^ ) : n £ No} are orthonormal basis that are dense in Ha 



and Hb, respectively. The sequence of states {<Jk}% x L\ in fact 
converges to oo m trace-norm, i.e., lim^oo || 0^ — colli = 
while liirifc^oo E(ak) = 00. Obviously, E is not continu- 
ous around the state <jq. In fact, the states with infinite von- 
Neumann entropy is trace-norm dense in state space |49J]. As 
can readily be verified, the set of pure states with infinite en- 
tropy of entanglement is also dense in the set of all pure states, 
and hence, the entropy of entanglement is almost everywhere 
infinite. That is to say, the assignment of a value which rep- 
resents the degree of entanglement to an infinite-dimensional 
quantum system is not entirely unambiguous. 

This apparent counterintuitive pathology can nevertheless 
be tamed to a significant degree. The key observation here 
is concerned with the mean energy of the states. Let H = 
Ha <8> i+ 1® Hb be the Hamiltonian of the bi-partite quantum 
system. It may appear to some extent unusual in the field of 
quantum information theory to refer to the Hamiltonians of 
the actual physical systems carrying the quantum states. We 
do not require much from H at all, however, only that for any 
finite temperature T we have 



tile- "} < 00 



(9) 



for all j3 > where (3 = l/T. This is a very natural de- 
mand on the spectrum on H. It merely means that there can 
be no limiting points in the spectrum of H: What it physi- 
cally implies is that the Gibbs state, the state of the canonical 
ensemble, exists. Systems such as the photon number degree 
of freedom of field modes of light, have this property. What 
are the implications for the discontinuity of the entropy of 
entanglement? We observe with Ha = Y1V=o^\^a)^a\ 

and Hb = Y^h=o M^B /i&B I tnat m tne arjove example the 
mean energy of the sequence of states diverges, 



lim ti[akH] = lim > 

t" — i rv-i t" — v rvi ■ " 



{ log(n) 2 



00. 



(10) 



This divergence of the mean energy is actually generic for se- 
quences of states with divergent entropy of entanglement 144], 
and is hence not an accident. This observation suggests how 
to tame the unwieldy infinities. After all, the energy that can 
be invested in the preparation of a state by means of physical 
devices is in all instances limited. This does not mean that 
the resulting states have a finite-dimensional carrier. It only 
means that their mean energy must be bounded. In fact, it 
turns out that for the set of states with mean energy which is 
bounded from above by a positive number E max , 



{p £ S(H) : H[ P H] < E max } 7 



(11) 



the ent ropy of entanglement regains its (trace-norm) continu- 
ity 1 44, 49]. Here, S(TL) denotes the state space, i.e., the set of 
all density matrices, which is the set of all positive normalized 
trace-class operators. The set Eq. ( II 11 is technically speaking 
nowhere dense in state space as H is an unbounded operator, 
but it is a reasonable subset of the state space: it simply re- 
flects the practically natural requirement that the mean energy 
is bounded from above. 
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Moreover, restricting the mean energy of the states un- 
der consideration recovers a number of asymptotic continu- 
ity properties of the entropy of entanglement. This strength- 
ens the interpretation as the entanglement cost and the dis- 
tillable entanglement also in the infinite-dimensional setting 
Similar conclusions can be reached for the entangle- 
ment of formation and the relative entropy of entanglement 
ll44ll - In a nutshell, one may say that on subsets of state space 
corresponding to bounded mean energy, entanglement mea- 
sures often regain trace-norm continuity. Nevertheless, one 
has to keep in mind that the actual minimization problems 
involved in the evaluation of meaningful entanglement mea- 
sures are often not feasible in the infinite-dimensional setting. 
This may again be taken as a very discouraging observation. 
However, in many situations of practical interest, where the 
semantics of the term practical may range from 'meaningful 
in mathematical physics' to 'preparable in a quantum optical 
experiment', it turns out that one does not encounter all pos- 
sible states from state space. Instead, a most relevant class of 
quantum states can be described in very simple terms using 
(small) finite-dimensional matrices only, without the need of 
a description that is overburdened with the technicalities of 
infinite-dimensional Hilbert spaces: this is the set of Gaussian 
quantum states. This set of states that is of utmost practical 
importance will be dealt with in most of the remainder of the 
present article. 



ENTANGLEMENT PROPERTIES OF GAUSSIAN STATES 

Gaussian quantum states play a key role in several fields of 
theoretical physics. In quantum optics, for example, they are 
often encountered as states of field modes of light, for reasons 
that will be elaborated on below. Ground states of systems 
with canonical coordinates (position and momentum) where 
the Hamiltonian is quadratic in the positions and momenta 
are also Gaussian states 1 50] . The term has been coined be- 
cause the defining property is that the characteristic function 
associated with the state is a Gaussian in phase space. For 
such Gaussian states the theory of entanglement as well as a 
framework of how these states may be manipulated is well de- 
veloped. Typical questions in a theory of entanglement, e.g., 
concerning the separability of given states, or the intercon- 
vertability of pairs of states under local operations can often 
be answered in full. This is true even in cases where the finite 
dimensional equivalent, whenever such an equivalent can be 
formulated, remains unsolved. One key reason for these suc- 
cesses is the fact that Gaussian states are completely specified 
by their first and second moments so that questions concerning 
properties of Gaussian states can be translated into properties 
of (comparatively small) finite-dimensional matrices. There- 
fore the instruments of matrix theory 12211 - so useful in the 
study of entanglement for finite-dimensional systems such as 
qubits and qutrits, once more become a useful tool. 

The systems that will subsequently be discussed are quan- 
tum systems with n canonical degrees of freedom. These 



could represent n harmonic oscillators, or n field modes of 
light. The canonical commutation relations (CCR) between 
the 2n canonical self-adjoint operators corresponding to posi- 
tion and momentum of such a system with n degrees of free- 
dom, may be written in a particularly convenient form em- 
ploying the row vector, 



— (Ol, ■ ■ . . C>2n) T — {X\. P\, ... , X n , P n 



(12) 



In terms of the familiar creation and annihilation operators of 
the modes choosing H = 1, X n and P n can be expressed as 

X n = (a n + 4)/V2, P n = -i(a n - 4)/\/2. (13) 

Then the canonical commutation relations (CCR) take the 
form 



[Oj,O k 



(14) 



where the skew-symmetric block diagonal real 2n x 2n-matrix 
a given by 



e 



1 

-1 



(15) 



is the so-called symplectic matrix. The phase space, isomor- 
phic to 1R 2 ™, then becomes what is known as a symplectic 
vector space, equipped with the scalar product corresponding 
to this symplectic matrix. Instead of referring to states, i.e., 
density operators, one may equivalently refer to functions that 
are defined on phase space. There are many common choices 
of such functions in phase space, such as the Wigner function, 
the Q-function or the P-function, to name just a few l5lll . 
Each of them is favorable in a particular physical context. For 
later purposes it is most convenient to introduce the charac- 
teristic function, which is the Fourier transform of the Wigner 
function. Using the Weyl operator 



(16) 



for £ e R , we define the (Wigner-) characteristic function 

as 



Xp(Z)=tr\pW 6 ]. 



(17) 



In quantum optics, the Weyl operator is typically referred 
to as phase space displacement operator or Glauber opera- 
tor, but with a different convention concerning its arguments. 
For a single mode, let the complex number a be defined as 
a = -(6 + *6)/V2, a* = -(£1 + then the phase 

space displacement operator D a of quantum optics |52] is 
most commonly taken as D a = W^. Each characteristic func- 
tion is uniquely associated with a state, and they are related 
with each other via a Fourier- Weyl relation. The state p can 
be obtained from its characteristic function according to 



1 



(270' 



(18) 



6 



In turn, the Wigner function as commonly used in quantum 
optics is related to the characteristic function via a Fourier 
transform, i.e., 



M0 = 



1 



(19) 



Gaussian states are, as mentioned before, defined through 
their property that the characteristic function is a Gaussian 
function in phase space [53], i.e., 



(20) 



where r is a 2n x 2n-matrix and D E R 2 ™ is a vector. As a 
consequence, a Gaussian characteristic function can be char- 
acterized via its first and second moments only, such that a 
Gaussian state of n modes requires only 2n 2 + n real param- 
eters for its full description, which is polynomial rather than 
exponential in n. The first moments form a vector, the dis- 
placement vector d £ R 2 ™, where 



d j = (O j ) p =tt[O j p], 



(21) 



j = 1 , . . . , 2n. They are the expectation values of the canonical 
coordinates, and are linked to the above D by D = ad. They 
can be made zero by means of a translation in phase space of 
individual oscillators. As a consequence the first moments do 
not carry any information about the entanglement properties 
of the state. The second moments are embodied in the real 
symmetric 2n x 2n covariance matrix 7 defined as 



7jlfc = 2Retr[p(O j - <0,) p ) (O k - (O 



k) P ) 



(22) 



With this convention, the covariance matrix of the n-mode 
vacuum is simply l2„. Again, the link to the above matrix T 
is r = cj T "f(j. Clearly, not any real symmetric 2n x 2n-matrix 
can be a legitimate covariance of a quantum state: states must 
respect the Heisenberg uncertainty relation. In terms of the 
second moments the latter can be phrased in compact form as 
the matrix inequality 



7 + ia > . 



(23) 



In turn, for any real symmetric matrix 7 satisfying the uncer- 
tainty principle d23l > there exists a Gaussian state the second 
moments of which are nothing but 7. So Eq. J23l > implies the 
only restriction on legitimate covariance matrices of Gaussian 
quantum states. 

This observation has quite significant implications concern- 
ing the question of separability of two-mode Gaussian states 
shared by two parties: it has earlier been pointed out that 
partial transposition is a positive, but not completely positive 
map. Hence, partial transposition must map separable states 
onto separable states, but there exist states for which the par- 
tial transpose is no longer positive. In fact, for bipartite qubit 
systems the positivity of the partial transpose is a necessary 
and sufficient criterion for separability. A necessary condition 
for separability of Gaussian states can be formulated imme- 
diately, once it is understood how partial transposition is re- 
flected on the level of covariance matrices. Indeed, it has been 



pointed out earlier that the partial transposition on qubits is in 
fact time-reversal. Time reversal in a system with canonical 
degrees of freedom is characterized by the transformation that 
leaves the positions invariant but reverses all momenta, 



X 



X Pi 



-P. 



(24) 



Let us now consider a system made up of 2n oscillators, where 
n are held by each party. Applying the time reversal operation 
to the n oscillators held by one of the parties, the covariance 
matrix will be transformed to a real symmetric matrix 7 given 
by 7 = FjF, with 



1 2 



e 



1 

-1 



(25) 



The matrix 7 is the matrix collecting the second moments of 
the partial transpose p TB of p. As the positivity of p Tn is 
equivalent to 7 satisfying the Heisenberg uncertainty princi- 
ple, one may now merely check whether 



7+ 4(7 > . 



(26) 



If 7 + ia is not positive, then the state p must in fact be en- 
tangled. For two-mode systems, one mode held by each party, 
the criterion ( I26> is both necessary and sufficient for sepa- 
rability of the Gaussian state. One direction of the proof has 
been exemplified above, the converse has been proven in Refs. 

Hence, the two-mode Gaussian state with covari- 
ance matrix 7 is separable exactly if the covariance matrix 7 
obtained via partial time-reversal in one sub-system also rep- 
resents a valid physical state, i.e., if 7 + ia > 0. 

It is remarkable that the above theorem can be extended 
considerably. In fact, for a bi-partite system the partial time- 
reversal also provides a necessary and sufficient criterion for 
the distillability of a bi-partite Gaussian state of two party 
holds one mode and the other party holds N modes, positiv- 
ity of the partial transpose is still a necessary and sufficient 
criterion for separability [l57ll . This result together with the 
discovery of bound entanglement in Gaussian systems rely on 
an alternative criterion for the separability of Gaussian states. 
In fact, a Gaussian state represented by the covariance matrix 
7 is separable if and only if there exist covariance matrices 7^ 
and 7b such that 



7 > 1A © 7b- 



(27) 



This formulation is a very useful tool to prove the validity of 
statements, but does not serve as a practical criterion for sep- 
arability, as it is generally not straightforward to find the ma- 
trices and 7s (but note that this may be checked in form 
of a semi-definite program). In general, given two parties, 
one holding M and the other N oscillators, both the question 
of separability 1 56, 58] and distillability can be decided effi- 
ciently, although the questions of distillability requires a more 
sophisticated approach |56]. 

These results are very encouraging but they come with one 
caveat. While the states have been restricted to be Gaussian, 
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no such restriction has been applied to the set of allowed op- 
erations. The result concerning the distillability for example 
only holds when arbitrary operations in infinite-dimensional 
systems are permitted, even practically very complicated ones 
such as the distinction between 999 and 1000 photons in field 
modes of light. In the light of this observation, one should 
not forget that much of the significance of the Gaussian states 
stems from the fact that the Gaussian operations are so impor- 
tant. This is by no means a tautology: The Gaussian opera- 
tions are those operations that change the state, but preserve 
the Gaussian character of any Gaussian input state. This set 
of operations is singled out from both a theoretical perspec- 
tive, but most significantly also from the perspective of prac- 
tical experimental implementation. A most relevant question 
is then what tasks can be accomplished with Gaussian states, 
but not under all physically available operations, but under 
Gaussian operations only. This observation motivates the ap- 
proach described in the next section, namely the investigation 
of what can and cannot be achieved when one restricts atten- 
tion to Gaussian states and the so-called Gaussian operations. 



CHARACTERIZING GAUSSIAN OPERATIONS 

Before we can begin the development of the theory of en- 
tanglement of Gaussian states under Gaussian operations, we 
need to characterize both Gaussian states (already done in the 
earlier sections) and Gaussian operations. The most straight- 
forward definition of Gaussian operations simply states that 
an operation is Gaussian exactly if it will map every Gaussian 
input state onto a Gaussian output state. The remainder of 
this section will present more practically useful characteriza- 
tions of Gaussian operations and in particular it will show that 
Gaussian operations correspond exactly to those operations 
that can be implemented by means of optical elements such 
as beam splitters, phase shifts and squeezers together with ho- 
modyne measurements I60LI61II - all operations that are in 
principle experimentally accessible with present technology 

Let us start with Gaussian unitary transformations. The 
most general real linear transformation S which implements 
the mapping S : O i — ► O' = SO will have to preserve the 
canonical commutation relations [Op 0' k ] 
the case exactly if the linear map S satisfies 



iojfcl. This is 



SaS 1 = 



(28) 



The set of real 2n x 2n matrices S obeying this condition form 
the so-called real symplectic group Sp(2n, M) l63l . whose el- 
ements are called symplectic or canonical transformations. To 
any symplectic transformation S also S T , S , —S are sym- 
plectic. The inverse of S is given by S^ 1 = aS T a^ 1 and 
the determinant of every symplectic matrix is dct[5] = 1 

El El. 

Every symplectic transformation corresponds to a 
unitary operation acting on the corresponding state space. As 
a consequence of the Stone-von Neumann theorem, given a 
real symplectic transformation S there exists a unique unitary 



transformation Us acting on the state space such that the Weyl 
operators satisfy UsW^Ug = W S £ for all £ e R 2 . On the 
level of covariance matrices 7 of an n-mode system a sym- 
plectic transformation S is reflected by a congruence 
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SjS 1 



(29) 



It is instructive to consider the generators of the Gaussian uni- 
tary transformations. These infinitesimal generators can be 
found by applying the condition that the canonical commuta- 
tion relations (on the level of Weyl operators) are preserved 
to an infinitesimal Gaussian unitary transformation. We con- 
sider U = e~ teG = 1 — ieG + o(e 2 ), for a real e and ar- 
rive at the conclusion that the generators G have the form 
G = J2jk=i9]k(OjO k + O k Oj)/2, which is a Hermitian 
quadratic expression in the canonical operators (and thus in 
the annihilation/creation operators). As a consequence, ev- 
ery Hamiltonian G can in quantum optical settings be imple- 
mented using beam-splitters, phase plates and mirrors only. 
Conversely, every array of passive optical elements is de- 
scribed by such a Hamiltonian. 

A particularly important subset of all symplectic transfor- 
mations is formed by those S £ Sp(2n, R) that are moreover 
orthogonal, i.e., K(n) = Sp(2n, R)nO(2ra). K{n) is again a 
group, and a compact one. The elements of K(n) correspond 
to the passive transformations |63|. In optical systems, these 
particular Gaussian unitary operations are those that can be 
realized using beam splitters and phase shifts, but excluding 
squeezers. Such a passive transformation does not change the 
eigenvalues of the covariance matrix and in an optical systems 
does not change the total photon number. Hence, it can never 
transform a state that is not squeezed into a squeezed one. A 
state is called squeezed, if one of the eigenvalues of the co- 
variance matrix 7 is smaller than unity 15311 . In more physical 
terms this means that using passive transformations only the 
state can be transformed into a state such that the uncertainty 
with position is smaller than the respective uncertainty of the 
vacuum state. Any S € Sp(2n, R) can be decomposed into 



S 



dj 
1/dj 



(30) 



with K, L 6 K(n) and dj E R. This is the Euler decompo- 
sition of symplectic transformations. Physically, this means 
that any Gaussian unitary transformation can be realized as a 
(i) passive transformation, a (ii) single-mode squeezing oper- 
ation on each of the n modes, and a (iii) subsequent second 
passive transformation. 

Noisy channels, such as those corresponding to optical 
fibers in optical systems, can only be modeled as an irre- 
versible quantum operation. A Gaussian channel 1 64, 
is such an irreversible quantum operation, p 1 — ► £(p). It is 
reflected by a transformation of the covariance matrix accord- 
ing to 



7 1 



AjA T + G. 



(31) 



8 



where A, G are 2n x 2n-matrices, and G is moreover sym- 
metric. The complete positivity of the quantum operation is 
reflected by the condition 



trix 



G + ia- iAaA T > 0. 



(32) 



Such channels emerge naturally if a system is coupled to de- 
grees of freedom that cannot be controlled. Often, this inter- 
action is of Gaussian nature, as the Hamiltonian is a polyno- 
mial of second degree in the canonical coordinates |67|]. As 
the environment degrees of freedom cannot be accessed, one 
considers a partial trace with respect to them. Optical fibers 
manifest themselves for example as Gaussian channels of this 
type for field modes |68]. Then, the overall dilation is re- 
flected by a Gaussian completely positive map with a trans- 
formation rule for the second moments according to Eq. Pill , 
Note that there are no restrictions on the matrix A: this means 
that if one aims to implement an (impossible) operation with 
a transformation of the form 7 1 — > AjA T , Eq. 13 H together 
with i32\ show the fundamental limitations imposed by the 
CCR on this operation. This is relevant in the context of, for 
example, quantum amplifiers, attenuators, or phase conjugat- 
ing mirrors approximately realising a time reversal operation. 
Also, quantum cloning 1 62, 17011 and secret sharing I7lll have 
been formulated as a Gaussian channel'. 

In turn, one may ask, whether it is true that any channel of 
this form can be realized in three steps: (i) appending an ad- 
ditional system prepared in a Gaussian state ('ancillae'), (ii) 
a Gaussian unitary operation including all modes, and (iii) a 
partial trace with respect to the environment ('discarding the 
ancillae')? The structure of such a dilation has been depicted 
in Fig. [2] Clearly, any such operation is reflected by a trans- 
formation of the form Eq. PH . The converse is not quite true, 
and one might have to 'add classical Gaussian noise', by dis- 
placing the state in phase space with a classical Gaussian dis- 
tribution. This is reflected by the addition of a positive matrix 
to the covariance matrix. 

Furthermore, one might after all perform selective measure- 
ments which preserve the Gaussian character of the state. Let 
us assume that n modes of light are entangled with a distin- 
guished n + 1st mode. If one projects the mode with label 
n + 1 onto the vacuum state with state vector |0), then the re- 
sulting state of the remaining n modes is still a Gaussian state. 
But what are its second moments? If the original covariance 
matrix is of the form 
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A C 
C T B 



(33) 



where A is an 2n x 2n-matrix and B is a 2 x 2 matrix cor- 
responding to the distinguished mode, then the resulting co- 
variance matrix can be determined by a Gaussian integration 



1 ^A-C{B + lL 2 „r i C J , (34) 
which is often referred to as the Schur complement of the ma- 



A C 
C T B + l 2 



(35) 



Note that this transformation is not of the form as specified in 
Eq. Pl> . Now, it follows how an ideal homodyne detection is 
reflected on the level of second moments. The resulting co- 
variance matrix does not depend on the outcome of the mea- 
surement |60|]. It is again a Schur complement of a matrix. 
Indeed, with it — diag(l, 0) we have 

7 1 — ► A — C(ttBit)~ 1 C t , 



(36) 



where the inverse is replaced by the Moore-Penrose pseudo- 
inverse, as the matrix is apparently no longer invertible. Eq. 
d36l gives the transformation law for the covariance matrix 
of n modes, after the mode with label n + 1 has undergone 
an ideal homodyne detection. In fact, any Gaussian operation 
can be thought of as a concatenation from ingredients of this 
above list l59l l60l l6lll . i.e., any physical operation mapping 
all Gaussians onto Gaussians can be thought of as a sequence 
of Gaussian channels (including Gaussian unitary operations), 
together with homodyne detection. This notion can be stated 
most clearly in terms of the isomorphism between positive 
operators and completely positive maps |72]. For Gaussian 
operations this isomorphism 1 59, 61] leads to the following 
form of a mapping of the general Gaussian operation, on the 
level of second moments, 



7^f 1 -f 12 (f 2 + 7 r 1 rT 2 , 

where we have denoted 

" ri r 12 

r 2 



r = 



1 12 



(37) 



(38) 



and r = FTF, with F being defined as in Eq. I25i with 
appropriate size. All Gaussian completely positive maps on 
n modes are characterized by a 2n-mode covariance matrix 
r > i<7 2 „ and a displacement vector d £ R 4 ™ 1 59] . Equipped 
with the transformation laws under these operations we may 
now turn to assessing what state transformations can be done 
in the Gaussian setting. 



MANIPULATION OF GAUSSIAN STATES UNDER 
GAUSSIAN OPERATIONS 

In the last section we have presented useful characteriza- 
tions of Gaussian operations. These characterizations form 
the essential tools for the study of the Gaussian local manip- 
ulation of Gaussian states and play important roles in many 
proofs in this area. We will now proceed with an exposition 
of known results in the area that have been discovered in the 
past few years. It should be noted however, that we will not 
proceed along the historical lines of development. While the 
program of the study of Gaussian states under Gaussian oper- 
ations has originated in Ref. 166] for mixed Gaussian states, 
we will begin with an exposition of the state of knowledge for 
pure states. 



9 



Pure state entanglement 

As can be expected, the most complete understanding of the 
behaviour of Gaussian states under Gaussian operations has 
been reached for pure states. Indeed, necessary and sufficient 
conditions for the possible state transformations under local 
Gaussian operations (GLOCC) can be given. Local means 
here that the applied Gaussian quantum operations are applied 
locally, accompanied by classical communication. This result 
crucially depends on a normal form that can be obtained for 
any bipartite system of continuous variables For 
any pure Gaussian state with covariance matrix 7 of an n x 
n-mode system, there exist local symplectic transformations 
Sa, Sb such that 



= © 



(S A © S B h(S A © S B ) T (39) 

~ cosh(2r fe ) sinh(2r fc ) 

cosh(2rfc) — sinh(2rfc) 

sinh(2r fc ) cosh(2r fc ) 

- sinh(2r fc ) cosh(2r fc ) 



k=l 



with rfc £ [0,oo). In other words, by means of local uni- 
tary Gaussian operations the pure Gaussian state described 
by 7 can be transformed into a tensor product of two-mode 
squeezed states 

00 

|Vfe) = (1 - (tanhr-fe/2) 2 ) 1 / 2 ^(tanhr fc /2)>, n) (40) 

71=0 

characterised by squeezing parameters r% (see Fig. As a 
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FIG. 1 : Schmidt decomposition of pure bi-partite Gaussian states. 

consequence the non-local properties of a pure Gaussian state 
of n x n modes can be characterised by a vector 



r = (n, 



(41) 



of two-mode squeezing parameters, which we will always as- 
sume to be given in descending order. This vector will play 
a role closely analogous to that of the Schmidt-coefficients in 
the majorization criterion for pure state entanglement trans- 
formations under general LOCC operations. The criterion for 
pure-state transformations under GLOCC can now be stated 
in a very simple manner in terms of this vector of two-mode 
squeezing parameters. We write r > r' iff ru > ri for all 
k = 1, . . . , n, in descending order. Then, the criterion for the 



possibility of transforming the initial to the final state can be 
expressed in an extraordinarily simple manner: we have that 



p' under GLOCC, iff r > r'. 



(42) 



Additional modes 




FIG. 2: Representation of a Gaussian channel. 

It should be noted that the criterion for pure state transfor- 
mations under Gaussian local operation has a simpler struc- 
ture than the criterion that determines which state transforma- 
tions are possible under general LOCC, not restricted to Gaus- 
sian transformations. As already stated in earlier sections, we 
have that 



P 



under LOCC iff X(ti A [p\) -< A(tr A [p']). (43) 



where -< denotes the majorisation relation Eq. (@). In other 
words, a pure state p can be transformed into another pure 
state p' if and only if the reduction to one part of the joint 
system is initially more mixed than finally in the sense of ma- 
jorisation theory. 

This criterion for finite-dimensional systems can be imme- 
diately carried over to the infinite-dimensional case, when ap- 
propriate care is taken when extending the definition of LOCC 
convertibility: one then writes a — ► p under LOCC if there 
exists a sequence of LOCC, such that the images of the maps 
approximate p in trace-norm. It is then instructive to com- 
pare the restrictions imposed by the majorisation criterion 
with the constraints under Gaussian LOCC. To start with, un- 
der LOCC, catalysis of entanglement manipulation is possible 
fTlil75ll : there exist pure states of finite-dimensional bi-partite 
quantum systems p, p' and ui such that 



P 



P 

) ui 



p' under LOCC, but 
p' © lj under LOCC . 



(44) 



The state ui serves as a "catalyst", as the entanglement inher- 
ent in the state is not consumed in the course of the trans- 
formation. Such an effect is not possible under GLOCC and 
pure Gaussian states, as collective operations are never more 
powerful than operations on pairs of quantum systems at a 
time. This notably refers to entanglement transformations 
only, as one can think of procedures estimating certain prop- 
erties of Gaussian states where collective operations are su- 
perior to product operations. It should moreover be noted 
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that the theorem governing interconvertability under GLOCC 
also covers the stochastic interconversion of pure Gaussian 
under GLOCC, as it has been pointed out in the previous sec- 
tion that any Gaussian transformation can be lifted to a trace- 
preserving operation. This is again in contrast to the finite- 
dimensional case 

It is also an immediate but remarkable consequence of 
the above criterion that under Gaussian local operations with 
classical communication, one cannot "concentrate" two-mode 
squeezing of pure Gaussian states. In particular, 



p{rf 



p[r') ®p(0)® (n_1) , underGLOCC, (45) 



for any n G N and any r' > r, where p(r) denotes a two-mode 
squeezed state with squeezing parameter r e [0, oo) . This is 
very much in contrast to the situation when general operations 
are allowed for: for any r £ [0, oo) there exists a r' > r such 
that 



p{rf 



p(r')®p(Q) under LOCC, 



(46) 



as has been shown in Ref. H73H . That is to say, this criterion 
already points towards the impossibility of distilling Gaussian 
states with Gaussian operations. 



Mixed states 

For mixed states, both the number of parameters character- 
izing the state increases and the structure of state transforma- 
tions becomes far more involved. To start with, a simple nor- 
mal form such as the Schmidt decomposition or the vector of 
two-mode squeezing parameters r is not available. The sim- 
plest normal form that can be obtained for two harmonic os- 
cillators is given by the Simon normal form for the covariance 
matrix 1 54] . If we are given a Gaussian state p of a two-mode 
system with covariance matrix 7, then there exist symplectic 
transformations Sa, Sb € Sp(2,M.) such that 
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(47) 



Further to the enlarged number of parameters characterizing 
the state, also the amount of classical communication in gen- 
eral Gaussian local operations cannot be bounded anymore as 
opposed to the pure state case where one round of one-way 
communication is sufficient. This latter fact has also been 
made use of in the proof of the key statement of Ref. Jl3l- 
Nevertheless, under certain restrictions on the available proto- 
cols and with the help of the Simon normal form it is possible 
to obtain some statements concerning the interconvertability 
of mixed states and one can find necessary and sufficient crite- 
ria for a transformation to be possible. One such restriction is 
that to local Gaussian channels as have been discussed above: 
they can be realized by adding an ancilla system, followed by 



a joint unitary transformation and the subsequent discarding 
of the ancilla (see Fig.|2]for illustration). 

For such local Gaussian channels, one may write 



P 



p' under LOG 



(48) 



for Gaussian two-mode states p and p', if there exists local 
Gaussian channels £a and £b such that 



(£a®£b)Q>)=P? 



(49) 



(LOG stands for local Gaussian channel). In fact, both neces- 
sary and sufficient conditions for the interconvertability under 
LOG can then be proven. For details we refer the reader to the 
literature, in particular Ref. |66] and for generalizations for 
example to a three-party version of the above statement |76]. 



Distillation of Gaussian states with Gaussian operations 

The availability of distillation protocols will be of crucial 
importance in the infinite-dimensional setting when we wish 
to implement long-distance quantum communication based on 
continuous variables. In fact, such distillation protocols must 
be performed at the beginning of any procedure that relies on 
the availability of highly entangled approximately pure shared 
entangled states. Noise due to unwanted coupling to an en- 
vironment is never entirely avoidable, and protocols have to 
be devised that effectively reverse this process. In the finite 
dimensional setting, distillation protocols are known, and us- 
ing the polarization degrees of freedom of light, they can be 
implemented optically (although they require photon counters 
in an iterative procedure that distinguish with great efficiency 
different numbers of photons) lf77ll . 

At the time it seemed fairly natural to expect that such dis- 
tillation schemes can also be constructed in the case of Gaus- 
sian states and Gaussian operations. Let us hence for a mo- 
ment assume that pure Gaussian two-mode squeezed states 
have been transmitted through lossy optical systems such as 
optical fibers. The resulting states are still Gaussian, provided 
that the noisy channel was a Gaussian channel, as is a good 
assumption in case of optical fibers. In order to distill the en- 
tanglement in an iterative procedure, one would take two pairs 
of identically prepared systems p®p and would feed them into 
one step of the procedure: the output could then be taken as 
the input of the next step. The states are identical if they un- 
derwent the same losses. The corresponding modes will from 
now on be denoted as Al, A2, Bl, and B2. A feasible itera- 
tive distillation protocol preserving the Gaussian character of 
an input state would look as follows (see Fig.[3j: 

(i) Application of a local Gaussian unitary operation, i.e., 
a map of the form 



p®p< — ► (U A <E> Ub){p®p){Ua®U b ) ] - 



(50) 



Note that the denotes the tensor product between dif- 
ferent copies, while the (g> denotes the tensor product 
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between the two parties. As has been discussed before, 
this set includes the passive transformations embody- 
ing those operations that can be implemented in optical 
systems using beam splitters and phase shifts. It how- 
ever also includes squeezing operations. For the two- 
mode case, the set of local unitary transformation is a 
20-dimensional manifold. Note that it is not required 
that both parties implement the same operation. 

(ii) A homodyne measurement on the modes A2 and B2. 

(iii) Finally, the parties communicate classically about the 
outcome of the measurement, and perform postprocess- 
ing of the states of modes Al and Bl with unitary Gaus- 
sian operations. 



for a state p, where || . || i again denotes the trace norm, and p TA 
is the partial transpose of p. For Gaussian states it has an inter- 
pretation in terms of a certain asymptotic entanglement cost 
II 611 . In the meantime, it should be noted, the entanglement 
of formation is also available for the case of symmetric two- 
mode Gaussian states, i.e., states for which the two reductions 
of the covariance matrix are identical up to local symplec- 
tic transformations 1 78 ] . For pure (and for symmetric mixed) 

Gaussian states the logarithmic negativity is related to the de- 

j i 

gree of squeezing in a monotone way (see, e.g., Ref. 1 79]). 
These considerations, hence, imply again that with Gaussian 
operations alone two identically prepared two-mode squeezed 
states cannot be transformed into a single two-mode squeezed 
state with a higher degree of squeezing, as has already been 
pointed out in Eq. i42i . 



Minimal non-Gaussian extensions for distillation 




FIG. 3: The class of considered feasible Gaussian distillation proto- 
cols. 

It was expected that Gaussian operations would be suffi- 
cient to provide for such an entanglement purification proto- 
col. It was however proven in Ref. 1 60] that this is not the 
case. Later this proof was extended to demonstrate that en- 
tanglement distillation of Gaussian state with Gaussian oper- 
ations is not possible in general, allowing in particular for a 
2n-mode input 1 59]. This is after all quite an astonishing re- 
sult: no matter how the Gaussian local operation with classical 
operation is chosen, the degree of entanglement can not be in- 
creased. Whatever operation is implemented, it will result in 
a loss rather than in an increase of entanglement. The opti- 
mal procedure is simple: not to even try to distill with Gaus- 
sian operations alone. This remarkable result also implies that 
quantum error correction is not feasible if one is restricted to 
Gaussian states and Gaussian operations. Hence, it becomes 
clear that is necessary to add some non-Gaussian resource to 
the repertoire to ensure that interesting quantum information 
processing tasks can be carried out in the presence of noise. 

The statement has been based on the quantification on the 
degree of entanglement in terms of the logarithmic negativity 
that has been mentioned before, which is defined as 

log 2 ||p T *||i (51) 



This difficulty sketched in this section - the impossibility of 
distilling Gaussian states with Gaussian operations - can how- 
ever be overcome. It has been known for quite some time that 
the unlimited availability of arbitrary non-Gaussian operations 
essentially requiring quantum computation allow for exam- 
ple for entanglement distillation of Gaussian states 145. 1.5611 to 
finite-dimensional singlets. This is an important result from 
the perspective of the theory of entanglement. Yet, not only 
from a practical point of view it would be crucially important 
to find solutions that make merely use of minimal extensions 
of the set of operations into the non-Gaussian regime, includ- 
ing only those additional operations that are to a high extent 
feasible in quantum optical settings. 

In a nutshell, it turns out that it is possible indeed to distill 
continuous-variable entanglement yielding highly entangled 
approximately Gaussian states, and this can be done in opti- 
cal systems using passive optical elements and photon detec- 
tors only. Such procedures have been introduced and investi- 
gated in detail in Refs. 1 80 , 81]. The key idea is to leave the 
Gaussian setting in a first non-Gaussian step which involves 
a single measurement only (corresponding to the positive out- 
come of a photon detector), and then apply Gaussian oper- 
ations only. The philosophy of the protocol is sketched in 
the following: Some source provides a two-mode squeezed 
state which is transmitted through fibers whose imperfections 
turn the states into symmetric noisy but still Gaussian states. 
At that stage a single non-Gaussian operation is carried out 
which takes the states out of the Gaussian state space and may 
perhaps also increase entanglement. Subsequent to that only 
Gaussian operations are applied in an iterative scheme which 
serve two aims. Firstly, they further increase the amount of en- 
tanglement and secondly, they make the states progressively 
more Gaussian. Asymptotically the states may then have a 
higher degree of entanglement than the original supply and at 
the same time approach Gaussian states as closely as possible. 

Surprisingly indeed, such a scheme exists using as the 
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FIG. 4: A possible scheme for entanglement distillation of Gaus- 
sian states using minimal non-Gaussian resources. Gaussian entan- 
gled states are corrupted during transmission but retain their Gaus- 
sian character. A single non-Gaussian operation is applied taking the 
states out of the Gaussian state space. Subsequently only Gaussian 
operations are applied which both distill entanglement and drive the 
states closer to the set of Gaussian states. 

only non-Gaussian resource a photo-detector that can distin- 
guish between the presence and absence of photons but is not 
able to count the precise number of photons (yes-no detec- 
tor) lliol lilll . The first non-Gaussian step takes two identi- 
cal copies of two-mode squeezed states, and measurements 
are performed on two of the four modes. The outcome is 
accepted as successful in case that both detectors click (i.e., 
register at least one photon), corresponding to the Kraus oper- 
ators E% = i — |0) (0|. The reflectivity and transmittivities of 
the beam splitter V are not 50 : 50, but have to be tuned ap- 
propriately. The following iterative Gaussian procedure takes 
two pairs p p of non-Gaussian states p as input which are 
then mixed at a 50 : 50 beam splitter, where the transforma- 
tion induced by the beam splitters are given by unitaries 

jj _ rpn 1 e ~R* a\ax e Ra 2 a\rp-n 2 (52) 

with T = R = l/\/2. Two of the output modes are then fed 
into a photon detector, each associated with Kraus operators 

E l = \0)(0\, £ 2 = 1-|0)<0|, (53) 

where |0) denotes the state vector associated with the vacuum 
state. The state is kept in case of the vacuum outcome of both 
local detectors. The unnormalized final state after one step is 
hence given by 

p = (Q,0\(U ®U){pQ p){U ® U) f \0,0). (54) 

The resulting two-mode states then form the basis of the next 
step. It is an iterative protocol, and it is event-ready, in the 
sense that one has a classical signal at hand which indicates 
whether the procedure was successful or not. No further post- 
processing has to be performed. Generic weak convergence 
to two-mode Gaussians can indeed be proven, and it can be 



demonstrated that the procedure often leads to highly entan- 
gled and squeezed states. Remarkably indeed, the scheme 
turns out to be fairly robust with respect to detector imper- 
fections 1 81]. It is beyond the scope of this paper to present 
details of this procedure, however, and the reader is referred 
to the literature here. 

Generation of entanglement with passive transformations 

The previous considerations were concerned with the lo- 
cal manipulation of Gaussian states. This is a meaningful ap- 
proach if one has a distributed system, where the necessity 
of local operations is simply dictated by the set-up. But how 
does one prepare the entangled states in the first place? Exper- 
imentally, there are several ways to do it, and again, it would 
be beyond the scope of this paper to present these possibili- 
ties (see, e.g., Ref. IIH). Instead, to exemplify the formalism 
developed before, we concentrate on the optimal creation of 
entanglement of Gaussian states using passive transformations 
only So in this subsection, the operations are 

allowed to be non-local with respect to the modes that have 
to be entangled. Passive operations are, again in optical sys- 
tems singled out as they correspond to those operations that 
can be realized with beam splitters and phase shifts. If simply 
the vacuum is fed into the additional input modes, then one 
even faces no practical problems of mode matching in optical 
systems. To start with, if one intends to entangle n Gaussian 
input modes with passive transformations, according to 

71 — ► 7' = S~fS T , S € K(2n), (55) 

the input states have to be squeezed. This is immediately ob- 
vious: any state that is entangled must be squeezed, since oth- 
erwise, the covariance matrix 7' of the n-mode output would 
satisfy 

7'>i2n©%n, (56) 

which in turn according to Hit implies separability. But what 
is the optimal operation, and, first and foremost, what is the 
optimal degree of entanglement that can be achieved in any 
two-mode output? 

This question can actually be completely solved, when the 
degree of entanglement is quantified in terms of the negativ- 
ity 17911 . Let 7 be the covariance matrix of the n-mode input 
state, Gaussian but not necessarily pure, and 7' = S^fS T be 
the resulting covariance matrix after application of the passive 
transformations. Then, the maximum amount of entanglement 
obtained for any two-mode sub-system is given by 

E N {p) = max (0, log(A 1 A 2 )/2) , (57) 

where Ai and A 2 are the two smallest eigenvalues of 7. This 
general solution of the problem establishes a generically valid 
link between the squeezing of the initial state to the degree 
of entanglement that can be potentially unlocked with pas- 
sive transformations. The proof of this statement as found in 
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Ref. ll79ll makes extensive use of the isomorphism between 
the groups K{n) and U(n). The optimal entangling trans- 
formation can moreover be constructively found, and in the 
practically important case of two input modes, a formula can 
be given for the appropriate choices of transmittivities, reflec- 
tivities, and phases of the optical elements. 

CONCLUSIONS 



ical appearance of classical properties in quantum Brownian 
motion |85], but also newly arising questions concerning the 
study of entanglement properties of interacting quantum sys- 
tems both in the static case 15011 and in their dynamical be- 
haviour 1 86, |87|. In these problems, the theory of Gaussian 
entangled states often allows for the exact analytical solution 
of many questions and therefore provides an ideal playground 
for the exploration of quantum entanglement. 



In this article we have briefly reviewed some elements of 
the theory of entanglement both in finite and infinite dimen- 
sional systems. We have formulated the basic questions that 
any theory of entanglement aims to answer when we consider 
entanglement as a resource for quantum information process- 
ing protocols. These key questions are those of characteriza- 
tion, (local) manipulation, and quantification of the entangle- 
ment resource. We then pointed out that entanglement theory 
in finite dimensional systems such as qubits generally requires 
the availability of arbitrary local operations. Statements con- 
cerning the feasibility or efficiency of a particular local state 
transformation may require the use of operations that are ex- 
traordinarily difficult to realize in practice. This motivated the 
development of a theory of entanglement under experimen- 
tally available operations and a particularly relevant example 
of such a theory can be found in infinite dimensional systems 
that are equivalent to harmonic oscillators. Examples of such 
systems are the photonic degree of freedom of a light mode, 
nano-mechanical oscillators and, approximately, cold atomic 
gases. The restricted set of experimentally available opera- 
tions that is being considered in this theory is that of Gaussian 
operations which allow for the creation and manipulation of 
Gaussian states. This restricted set of operations is of par- 
ticular importance in quantum optics where it can be shown 
that Gaussian operations are in one-to-one correspondence to 
the set of operations that can be implemented by simple op- 
tical tools such as phase plates, beam-splitters, and squeezer 
together with the addition of vacuum modes and homodyne 
detection. In the remainder of the article we then outlined 
the recent development of this theory and stated some of its 
key results. In particular it became clear that processes such 
as quantum error correction, entanglement distillation and ef- 
ficient quantum computation are not possible under this re- 
stricted set of operations but require some resources that lie 
outside the Gaussian regime. This suggested to search for 
protocols that employ non-Gaussian resources in a minimal 
way. Indeed, such extensions are possible in the sense that the 
expensive non-Gaussian operations are used only in an initial 
step of the protocol and are then followed by purely Gaussian 
operations. 

Finally it should be mentioned that the recent development 
of the theory of entanglement for Gaussian states has also 
provided us with many novel and useful techniques for the 
analytical study of entanglement in harmonic systems. This 
permits both, to revisit old questions such as for example con- 
cerning the connection between entanglement and the dynam- 
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